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Abstract 



The e + p neutral-current deep inelastic scattering differential cross-sections da/dQ 2 , for 
Q 2 > 400 GeV 2 , dcr/dx and da/dy, for Q 2 > 400, 2500 and 10000 GeV 2 , have been measured 
with the ZEUS detector at HERA. The data sample of 47.7 pb _1 was collected at a center-of- 
mass energy of 300 GeV. The cross-section, da/dQ 2 , falls by six orders of magnitude between 
Q 2 = 400 and 40000 GeV 2 . The predictions of the Standard Model are in very good agree- 
ment with the data. Complementing the observations of time-like Z° contributions to fermion- 
antifermion annihilation, the data provide direct evidence for the presence of Z° exchange in 
the space-like region explored by deep inelastic scattering. 
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1 Introduction 



The HERA ep collider has made possible the exploration of deep inelastic scattering (DIS) in a new 
kinematic region, resulting in discoveries such as the rapid rise of the parton densities in the proton 
at low x-Bjorken Q and diffraction in DIS First measurements of charged and neutral-current 
(NC) deep inelastic electron-proton scattering were made in a previously unexplored region of large 
Q 2 Q, where Q 2 is the negative square of the four-momentum transfer between the electron and 
proton. Based on their 1994-1996 data, both the HI and ZEUS collaborations have reported [||, [5|] 
more NC events than expected from the Standard Model at high Q 2 and high x. These observations 
have prompted considerable discussion in the particle-physics community as possible evidence for 
anomalies in the parton momentum distributions of the proton or of physics beyond the Standard 
Model. 

This paper presents measurements at HERA of the NC deep inelastic positron-proton scattering 
differential cross-sections da/dQ 2 for Q 2 > 400 GeV 2 and da/dx and dajdy for Q 2 > 400, 2500, and 
10000 GeV 2 , and their comparison to Standard Model predictions. This analysis is not optimized 
for the search for narrow lepton-hadron resonances. The measurements are based on 47.7 pb _1 
of data collected by ZEUS from 1994-1997 during which HERA collided 27.5 GeV positrons with 
820 GeV protons, yielding a center-of-mass energy yfs = 300 GeV. The highest Q 2 under study, 
51200 GeV 2 , is much larger than the square of the Z°-boson mass so that effects of Z° exchange 
are visible. 



2 Standard Model prediction 

The electroweak Born-level NC DIS differential cross-section, d 2 a^ in /dxdQ 2 , for the reaction 
e+p -> e + X can be expressed || as 

d2a 2 r dQ 2+P) = i$ [Y+ F * C{X > Q2) ~ Y - xF ' C{x > Q2) ~ y2 F " C{X > Q2)] ' (1) 

where x is the Bjorken scaling variable, a(Q 2 = 0) ~ 1/137 is the QED coupling constant, and 
Y± = 1 i (1 — y) 2 with y = Q 2 /sx. The structure functions F^ c and xF^ c for longitudinally 
unpolarized beams may be described in leading order QCD as sums over the quark flavor f = u, ...,b 
of the product of electroweak quark couplings and quark momentum distributions in the proton 

F 2 NC = I £ xg+ [(U/) 2 + (Vf) 2 + (Ajf + (Af) 2 ] , 
/ 

^3 NC = E *9J [Vf A f ~ V f RA f] > ^ 
/ 

where xq^ = xqj(x,Q 2 ) ± xq~j(x,Q 2 ) and xqf (xqf) are the quark (anti-quark) momentum dis- 
tributions. In leading order QCD, we have F^ c = 0. The functions Vf and Af can be written 

as 

V f ,R = e f ~ ( V e ± a e) Vf Xz(Q 2 ), 

( 3 ) 

Af R = -(v e ±a e )a f Xz (Q 2 ), 
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where the weak couplings, aj = Tf and Uj = Tf — 2ejsin 2 6' w , are functions of the weak isospin, 
Tf = \ (—5) for u, z/ (d, e), and the weak mixing angle, 6 W ; e« is the electric charge in units of the 
positron charge; and x z i s proportional to the ratio of Z°-boson and photon propagators 

x = I 9L_ (4) 

Xz 4sin 2 # w cos 2 # w Q 2 + M 2 " {> 

All cross-section predictions in this paper are calculated using next-to-leading order (NLO) QCD 
where the longitudinal structure function F^ c 7^ @. The contribution of F^ c to d 2 a^ rn /dxdQ 2 
is predicted to be approximately 1.5% averaged over the kinematic range considered in this paper. 
However, in the region of small x at the lower end of the Q 2 range the F^ c contribution to the 
cross-sections can be as large as 12%. 

Uncertainties in the predicted cross-section arise from three sources: electroweak parameters, elec- 
troweak radiative corrections, and the parton momentum distributions including their higher order 
QCD corrections. The electroweak parameters have been measured to high precision by other ex- 
periments & P] and contribute less than 0.3% uncertainty to the predicted cross-section in the full 



kinematic range measured at HERA [1C]. Radiative corrections for initial- and final-state radiation, 



vertex and propagator corrections, and two-boson exchange have been calculated to at least first 



order [11, 12 1. Higher order corrections for the kinematic region explored in this paper are expected 



to be less than 1% [13]. This leaves the parton momentum distributions as the primary source of 



uncertainty in the predicted cross-section. 



Parton momentum distributions have been determined by several groups (GRV [14], MRS |L^, 16], 
CTEQ flT] , |18|| ) by parameterizing the distributions at some fixed Q 2 and extrapolating the results 
to higher Q 2 using the DGLAP QCD evolution equations The parameters are fitted to 

data from lower energy fixed-target DIS experiments and from HERA, and, in addition, to data 
measured at the TEVATRON on lepton-pair production (Drell-Yan), direct photon production, W 
production, and jet cross-sections. Note that the HERA data included in these parameterizations 
make their most significant contribution at x < 10~ 2 and have relatively little influence on the 
predicted cross-sections used in this paper. The sources of uncertainty in these fits can be divided 
into two main groups: uncertainties in the measurements and uncertainties in the fit itself. For 
the former, the statistical and systematic uncertainties are available from each experiment. For the 
latter, uncertainties due to the QCD evolution were estimated by varying the fit assumptions, such 
as the value of a s and higher twist. 

Although the comprehensive parton momentum distribution fits from GRV, MRS, and CTEQ make 
extensive use of available data, they lack a complete estimate of uncertainties in the distributions. 



To obtain such an estimate, a NLO QCD fit was performed [20] to the DIS measurements of 



F 2 for proton and deuteron data from SLAC jgy, BCDMS @, NMC f23|, E665 j2|], HI and 



ZEUS [^, ^6[, the measurements of xF% from neutrino measurements by CCFR [|27j, and the d — u 
data from E866 [28]. Included in the fit were statistical and correlated systematic errors from 



each experiment. Also considered were the effects of a change of a s (M|) from 0.113 to 0.123, 
a 50% variation in the strange quark content of the proton, a variation of the factorization and 
renormalization scales /j,f >r in the range Q 2 /2 < fi 2 r < 2Q 2 , and of corrections for nuclear effects, 
all of which produced uncertainties of typically less than 1%. The results are in good agreement 
with the MRST |l(| and CTEQ4 |l8| fits, and differences are typically smaller than the extracted 
uncertainties. The fit yields uncertainties in the cross-section da/dQ 2 of approximately 2.5% for 
Q 2 = 400 GeV 2 and 6% at the highest Q 2 under study. 
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Other uncertainties were also investigated. For example, charm-threshold effects, calculated from 
three different models of charm evolution in the proton as supplied by CTEQ |29f| , produced cross- 
sections that differed by less than 3%. An analysis of the stability of perturbative calculations of the 
production of bottom- and charm-quarks at HERA |3(]] showed negligible effects. The MRST [16] 
fit incorporating the latest prompt photon data at high x from E706 |H]] and d — u data from 
E866 p§ produced cross-sections that were lower than those calculated using CTEQ4 by 4% at 
Q 2 = 400 GeV 2 , increasing to a maximum of 8% at Q 2 = 10000 GeV 2 . The CTEQ4HJ @ fit, 



specifically tuned to reproduce the jet high transverse energy cross-section reported by CDF [32 



produced changes in the cross-section of less than 2% except at the highest Q 2 = 50000 GeV 2 , where 
it yielded an increase of 6%. The CTEQ5 fit incorporating more data than the CTEQ4 fit, in 
particular introducing an improved handling of d/u and d — u using E866 data and a measurement 
of charge asymmetry in W-production at the TEVATRON [34], gave cross-sections that were higher 
by 3% at Q 2 = 400 GeV 2 and lower by 2% at Q 2 = 10000 GeV 2 . 

We conclude from these studies that the parton densities give a total uncertainty on the Standard 
Model prediction of the NC DIS differential cross-section da/dQ 2 of 4% for Q 2 = 400 GeV 2 in- 
creasing to 8% at the highest Q 2 in the x-range covered by this measurement. In the following, 
uncertainties in the parton momentum distributions are taken from the ZEUS NLO QCD fit [20]. 



3 The ZEUS experiment 

ZEUS |5E| is a multipurpose magnetic detector designed to measure ep interactions at HERA. The 
primary components used for this analysis are the compensating uranium-scintillator calorimeter 
(CAL), the central tracking detector (CTD), and the luminosity detector. 

The ZEUS coordinate system is right-handed with the Z axis pointing in the direction of the proton 
beam (forward) and the X axis pointing horizontally toward the center of HERA. The polar angle 
9 is zero in the Z direction. 

The CAL covers 99.7% of the total solid angle. It is divided into three parts with a correspond- 
ing division in 9 as viewed from the nominal interaction point: forward (FCAL, 2.6° < 9 < 36.7°), 
barrel (BCAL, 36.7° < 9 < 129.1°), and rear (RCAL, 129.1° < 9 < 176.2°). Each section is sub- 
divided into towers which subtend solid angles between 0.006 and 0.04 steradian. Each tower is 
longitudinally segmented into an electromagnetic (EMC) and one (RCAL) or two (FCAL, BCAL) 
hadronic sections (HAC). The electromagnetic section of each tower is further subdivided trans- 
versely into two (RCAL) or four (BCAL, FCAL) cells. Under test beam conditions the calorimeter 
resolutions are a/E = 18%/-^/ E (GeV) for electrons and a/E = 35%/ y 7 E (GeV) for hadrons. The 
calorimeter has a time resolution of better than 1 ns for energy deposits above 4.5 GeV. 

A presampler detector is mounted in front of FCAL and RCAL. It consists of scintillator tiles 
matching the calorimeter towers and measures signals from particle showers in the material between 
the interaction point and the calorimeter. 

Tracking information is provided by the CTD |37j operating in a 1.43 T solenoidal magnetic field. 
The interaction vertex is measured with a typical resolution along (transverse to) the beam direc- 
tion of 0.4 (0.1) cm. The CTD is used to reconstruct the momenta of tracks in the polar angle 
region 15° < 9 < 164°. The transverse momentum (pt) resolution for full-length tracks can be 
parameterized as a{pt)/pt = 0.0058 p t © 0.0065 © 0.0014/p t , with p t in GeV. 
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The luminosity is measured using the Bethe-Heitler reaction ep — ► eyry p8| . The resulting small 
angle energetic photons are measured by the luminosity monitor, a lead-scintillator calorimeter 
placed in the HERA tunnel 107 m from the interaction point in the positron beam direction. 



4 Monte Carlo simulation 

Monte Carlo simulations (MC) are used to determine the efficiency for selecting events, to determine 
the accuracy of kinematic reconstruction, to estimate the background rate, and to extrapolate the 
measured cross-sections to the full kinematic phase space. A sufficient number of events is generated 
to ensure that errors from MC statistics can be neglected. The MC samples are normalized to the 
total integrated luminosity of the data. 

The ZEUS detector response is simulated with a program based on GEANT |3i|. The generated 
events are passed through the simulated detector, subjected to the same trigger requirements as 
the data, and processed by the same reconstruction programs. 

The vertex distribution is a crucial input to the MC as this is necessary to estimate the event 
selection efficiency. The latter is strongly correlated with the Z-coordinate of the event vertex, as 
the reconstructed values of Q 2 , Et and other quantities depend on the measured vertex position. 
For the 1995 to 1997 data sets, the underlying distribution of the Z-coordinate of the event vertex is 
determined using a minimum-bias sample of low Q 2 neutral-current DIS events. For 1994 data this 
method is compared to the estimate from a special minimum-bias soft photoproduction trigger, 
where very good agreement in shape is found. The uncertainty in the shape of the vertex Z- 
distribution is related to the fraction of ep collisions from RF buckets adjacent to that containing 
the main proton beam. The effect of this uncertainty on the normalization of the data sample in 
the vertex range -50 to +50 cm is estimated to be less than 0.4%. 



NC DIS events including radiative effects are simulated using the HERACLES 4.5.2 [11] program 



with the DJANG06 2.4 [ 40 1 interface to the hadronization programs. In HERACLES, corrections 
for initial- and final-state radiation, vertex and propagator corrections, and two-boson exchange 
are included. The QCD cascade and the hadronic final state are simulated using the color-dipole 



model of ARIADNE 4.08 pj and, as a systematic check, with the MEPS model of LEPTO 6.5 {42 



Both programs use the Lund string model of JETSET 7.4 p3] for the hadronization. 



Photoproduction background is estimated using events simulated with HERWIG [44]. In addition, 
a large sample of prompt photon events (ep — ► e'jX), is generated with HERWIG. 



5 Reconstruction of kinematic variables 

5.1 Event characteristics 

Neutral-current DIS at the high-Q 2 values discussed here produces striking events, relatively easy 
to distinguish from the potentially large backgrounds of quasi-real photoproduction (Q 2 ~ 0) and 
beam-gas interactions. The events are characterized by a high-energy isolated positron in the 
detector. For Q 2 > 400 GeV 2 , most of the positrons have an energy near the positron beam energy 
and are restricted to a polar angle below 140°. As Q 2 increases, the positrons are produced with 
higher energies, up to several hundred GeV, and at smaller polar angles. 
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The variables 5, Pt and Et are used for event selection. These are defined as 

5 = ^(^-^cos^) = ^(E-^)i, (5) 

i i 

where the sum runs over all calorimeter energy deposits E^ (uncorrected in the trigger, but corrected 
energies in the offline analysis as discussed below) with polar angle 0j. At the generator level 5 = 
55 GeV, i.e. twice the positron beam energy, which follows from energy-momentum conservation. 
Undetected particles which escape through the forward beam hole give a negligible change in 
5 while particle loss through the the rear beam hole, e.g. from initial state bremsstrahlung or 
for photoproduction background, can lead to a substantial reduction of 5. The net transverse 
momentum, Pt, and the transverse energy, Et, are defined by 




E T =Y,Ei sin U (7) 

i 

where fa is the azimuthal angle and the sums run, as above, over all energy deposits in the calorime- 
ter. High-Q 2 events in which the positron strikes the BCAL or FCAL (Q 2 > 1000 GeV 2 ) are 
characterized by large Et- 

In the determination of the DIS kinematics, the CAL energy deposits are separated into those 
associated with the identified scattered positron, and all other energy deposits. The sum of the 
latter is referred to as the hadronic energy. 



5.2 The double-angle method 

Q 2 , x, and y are measured using the double-angle method |45| 

n2 1F 2 sin 7 ft(l + cos0 e ) 

sin 7/! + sm0 e - sm( 7 / l + e ) 
E e sin7 fe + sin0 e + sin( 7fe + e ) 
DA E p sin j h + sin 9 e - sin( 7/l + 9 e ) ' 
sin0 e (l - cos 7 ft ) 

1/dA = -: —^-5 . "n~\1 (10) 

sm7ft + sm0 e - sm( 7/l + e ) 

where E e (E p ) is the energy of the positron (proton) beam, e is the polar angle of the scattered 
positron, and 7 ft, in the naive quark parton model picture of DIS, is the polar angle of the struck 
quark. The determination of the angles e and 7 ft is discussed in Sect. |T3] and ^J, respectively. 

The double-angle method is essentially insensitive to errors in the absolute energy scale of the 
calorimeter. However it is sensitive to QED radiation and an accurate simulation of the hadronic 
final state is necessary. At Q 2 > 400 GeV 2 the hadronic final state is sufficiently energetic that 
model uncertainties of fragmentation and the treatment of the proton remnant are less important 
than at lower Q 2 . 

To validate the performance of the double-angle method, the reconstructed kinematic variables of 
MC events are compared to the true hadron variables Q 2 , x and y as defined by the four- momentum 
transfer q to the hadronic system 

Q 2 = -q\ x = Q 2 /(2p-q), y = Q 2 /{xs) (11) 
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where p is the four-momentum of the initial proton and s = AE p E e . The resolution in the kinematic 
variables is determined accordingly (see Sect. 0) and demonstrates (not shown here) that the double- 
angle method performs better than other methods for Q 2 > 400 GeV 2 . 



5.3 Positron identification 

A key signature of the events under study is the presence of an isolated high-energy positron. In 
order to identify and reconstruct the positron, while rejecting events in which other final state 
particles mimic it, an algorithm is used which combines calorimeter and CTD information. 

In a first step, calorimeter "clusters" are formed by grouping all CAL cells with energy depositions 
into cones relative to the event interaction point, which are centered around cells with a local energy 
maximum [|^|. All clusters are treated as positron candidates. The cluster energy, E c \ u , is the sum 
of the cell energies belonging to the cluster. The center of each cluster is calculated by weighting 
each tower member by a logarithmic function of its energy. The cluster angle, $ciu? is set equal 
to the polar angle obtained from the center position of the cluster and the event vertex obtained 
from the tracks measured with the CTD. For clusters with polar angleQ within the CTD acceptance 
(#ciu > 17.2°), a matching track is required. A track is considered to match if the distance of closest 
approach (DCA) between the extrapolation of the track into the calorimeter and the position of 
the cluster center is less than 10 cm, where the r.m.s. resolution in the DCA is 1.8 cm. 

In the second step, several quantities, £j, are calculated for each positron candidate: the fraction 
of the cluster energy in the hadronic sections of the calorimeter; the parameters related to lateral 
energy profiles; and the total energy in all calorimeter cells not associated with the cluster but 
lying within a cone in pseudorapidity, r] = — ln(tan(0/2)), and azimuthal angle, $, of radius 
R = yj (Ar/) 2 + (A<£) 2 = 0.8, centered on the cluster. If a matching track is present, the polar and 
azimuthal angle differences between the track and the cluster position, and the quantity 1/E C \ U — 
1/Ptrki where P t rk is the track momentum, are evaluated. 

Finally, for each £j a probability P(£i) is derived, designed to be uniformly distributed between zero 
and one for positrons. Candidates are accepted as positrons if the product of the P(£i) exceeds a 
threshold determined from MC studies. Since the CAL energy resolution is better than that of the 
CTD for tracks above 10 GeV, for accepted candidates, the positron energy, E' e , is set equal to the 
cluster energy, E c \ n , comprising typically six cells. The positron angle, 6 e , is determined from the 
associated track if the positron cluster is within the CTD acceptance, and otherwise set equal to 
9 c i n . The resolution in 9 e is typically 0.1°. Sect. 6^2 contains further selection criteria applied to 
the positron candidates. 



5.4 Reconstruction of the hadronic final state 

Energy depositions of the hadronic final state are used to evaluate the angle jh- F° r this purpose 
calorimeter clusters are used as just described. The angle jh is then calculated for the event 

x We do not consider candidates with 6 c i u > 164° (which are also beyond the CTD acceptance limit), since they 
correspond to Q 2 values below the range of this analysis. 
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according to 



T,h °h nr> s 

cos 7 fe = 2 , 2 , (12) 
T,h °h 

where Pr,h an d $h are calculated using (||) and (||) with sums running over the calorimeter clusters 
in the same manner as Pt and 5, but excluding the positron. 



MC studies of the calorimeter response indicate that the uncorrected jh calculated with (12) is 
biased by redirected hadronic energy from interactions in material between the primary vertex and 
the calorimeter or by backsplash from the calorimeter (albedo )0. To minimize this bias, clusters 
with energy below 3 GeV and with polar angles larger than 7 max are removed. The value of 7 maX) 
which is a function of 7^, is derived from a NC MC sample by requiring that less than 1% of the 
clusters not related to the above effects be removed. This yields a reconstruction of 7^ closest to the 
true value as given by the MC. The agreement of the distribution of removed energies for different 
7/j-ranges between data and MC is reasonable as shown in Fig. |l](a) and (b). After this first pass 
of cluster removal the value of 7^ is re-calculated and the procedure is repeated until it converges, 
typically after two or three passes. Removing calorimeter clusters in this manner substantially 
improves the resolution and bias of the double-angle variables for small values of 7^ (corresponding 
to small values of y) as shown from the ARIADNE MC in Fig. |l](c) and (d) and leaves them largely 
unchanged for large values of jh- A similar result is achieved with the LEPTO MC (not shown 
here). 



5.5 Energy correction and calorimeter performance 

Correction for energy loss: All energies of clusters from both the positron and the hadronic 
final state are corrected for energy loss in the material between the interaction point and the 
calorimeter. All hadronic clusters and the positron clusters in BCAL and FCAL are corrected 
based on the material maps implemented in the detector simulation package. The presampler is 
used to correct positrons entering the RCAL. This correction is based on the measured positron 
energy, the amount of material in front of the RCAL and the presampler signal. 

Non-uniformity correction for positrons: In Fig. ^(a) the ratio of the positron energy 
corrected as described above to the double-angle prediction is shown. At boundaries of calorimeter 
cells and modules there are dips in the measured energy response. This effect is larger in data than 
in MC. To account for this effect a correction is obtained from fits to the non- uniformity patterns 
in the BCAL. After correction, the data are well reproduced by the MC as shown in Fig. ||(b). A 
similar correction is used for geometrical effects in the azimuthal angle. 

Calorimeter energy scale: The uncertainty in the energy scale of the scattered positron is 
determined after applying the corrections described above. For scattered positrons in low y DIS 
events [9 e > 135° due to the requirement y < 0.04), the energy is strongly correlated with the 
scattering angle, and a comparison of the predicted energy to the measured energy in the calorimeter 



is made. This method is described in detail elsewhere [26|. In the range of 30° < 8 e < 150°, the 



momentum of the positrons can be determined by the CTD. The average track momentum minus 



2 High energy hadrons interacting in a calorimeter can produce with a non- negligible probability particles at large 
angles with respect to the direction of the main shower. Some of these particles travel backwards and generate energy 
deposits far away from their primary source in the calorimeter which is referred to as backsplash. 
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calorimeter energy is used as an independent check of the energy measured in the CAL for energies 
up to 30 GeV. For positrons with energies above 30 GeV, or those scattered to extreme forward 
directions, a comparison of the energy predicted from double- angle variables to the measured energy 
is made. In kinematic regions where the other methods can be used to check the CAL energy scale, 
the double-angle results are in agreement with the other methods to better than 0.5%. As a result 
of these studies, the uncertainty in the value of the energy of the scattered positron in the RCAL 
is 2% at 10 GeV, decreasing linearly to 1% at 27.5 GeV and above, 1% in the BCAL and 3% in the 
FCAL. 

The uncertainty in the scale of the hadronic energy has been studied. After applying the correc- 
tions to the energy described above, the ratio of the hadronic transverse momentum, Pt,Hi to the 
transverse momentum Pr,e carried by the positron found in the BCAL, is examined event by event. 
For transverse momenta above 20 GeV, the ratio is 1.0 to within 0.5% both for data and MC, 
as expected from momentum conservation (see Fig. |3|). For lower transverse energies, the ratio is 
below 1.0 by up to several percent due to hadronic energy loss. The comparisons between MC 
predictions of the transverse momentum ratio and the data agree to within 1-2%. A comparison 
with corresponding results using the ratio of Pt,h to the transverse momentum predicted from 
double-angle variables also shows good agreement. As a result of these studies, the uncertainty in 
the hadronic energy scale is determined to be 2% in the FCAL and BCAL, of which 1% comes 
from the uncertainty of the positron measurement. The dominant remaining uncertainty comes 
from comparison of MC results for simulations of hadronic final states using ARIADNE |4l| and 
HERWIG [Q. The uncertainty is 3% in the RCAL where the hadronic energy is lower than in 
FCAL and BCAL due to kinematic constraints. 

Positron energy resolution: After all corrections, the resolution in positron energy is evaluated 
by comparing the width of the distribution of the ratio of corrected energy to the value predicted 
from the double-angle method. Since the resolution is smaller in the MC than in data, the MC 
energies in BCAL and RCAL are smeared accordingly. For RCAL a constant relative smearing of 
1.7% is applied while for BCAL the dependence of the resolution on the calorimeter non-uniformity 
is taken into account. 

Hadronic energy resolution: A similar method, based on Ptk-, is used to determine the 
resolution in hadronic energy. Again, the resolution is found to be smaller in the MC than in data. 
A smearing is applied accordingly to clusters in MC for all calorimeter parts. 

5.6 Detector alignment 

The polar and azimuthal angles of the scattered positron can be measured with the tracking detec- 
tors as well as with the calorimeter. Comparisons of the position of calorimeter non-uniformities 
resulting from the cell structure (see Sect. f>.5\) and the extrapolated positron position from the 
CTD allow the BCAL and RCAL to be aligned in Z with respect to the CTD to 1 and to 2 mm, 
respectively, and to be aligned in transverse direction to 2 and to 1 mm, respectively. 

6 Event selection 
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6.1 Trigger 

The ZEUS trigger operates at three levels |3j|. For the portions of the trigger relevant to this 
analysis, the requirements were strictest during 1997 and are described here. The first-level trigger 
requires a total electromagnetic energy of at least 3.4 GeV in the RCAL or 4.8 GeV in the BCAL, 
or E'Jp > 30 GeV, where Ell is the total transverse energy excluding the two rings of calorimeter 
towers nearest to the forward beampipe. The Ej, requirement is designed to tag high-Q 2 events 
by their large Ej~ while rejecting beam-gas background having large energy deposits at low polar 
angles. The major requirement at the second-level trigger is 5 + 2E^ > 29 GeV, where J5~ is 
the energy measured in the luminosity monitor. This requirement suppresses photoproduction. 
Backgrounds are further reduced at the second-level trigger by removing events with calorimeter 
timing inconsistent with an ep interaction. For the third-level trigger, events are reconstructed on a 
computer farm, and the requirements are similar to the offline cuts described below, though looser 
and using a simpler and generally more efficient (but less pure) positron finder. 

6.2 Offline selection 

The following criteria are applied offline (see also Sect. || and ]5[): 

• To ensure that event quantities can be accurately determined, a reconstructed vertex with 
—50 < Z < 50 cm is required, a range consistent with the ep interaction region. 

• To suppress photoproduction events where the scattered positron escapes through the beam 
hole in the RCAL, 5 is required to be greater than 38 GeV. This cut also reduces the number 
of events with initial-state QED radiation. The requirement 5 < 65 GeV removes cosmic ray 
background. 

• Positrons are identified based on calorimeter cluster quantities and tracking. 

— To ensure high purity, the positron is required to have an energy of at least 10 GeV; in 
this case the identification efficiency exceeds 96%, as shown by MC studies. 

— To reduce background, isolated positrons are selected by requiring no more than 5 GeV 
in calorimeter cells not associated with the scattered positron in an 77 — $ cone of radius 
0.8 centered on the positron. 

— In addition, each positron with 8 e > 17.2° must be matched to a charged track of at 
least 5 GeV momentum. 

— For positrons beyond the forward tracking acceptance (6 e < 17.2°), the tracking require- 
ment in the positron selection is replaced by a cut on the transverse momentum of the 
positron Pt e > 30 GeV and by the requirement 5 > 44 GeV. 

— A fiducial volume cut is applied to the positron position. This excludes the upper part of 
the central RCAL area (20 x 80 cm 2 ) occluded by the cryogenic supply for the solenoid 
magnet as well as the transition region between the CAL parts corresponding to a polar 
angle of the positron of 35.6° < 6 < 37.3° and 128.2° < 9 < 140.2°. 

If more than one positron candidate in an event passes these cuts (7% of the events), the one 
with the highest probability is assumed to be the DIS positron. 
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• Elastic Compton scattering events (ep — ► e^p) are removed by searching for an additional 
photon candidate and discarding the event if this and the positron candidate account for all 
the calorimeter energy in the event except for at most 2 GeV. 

• To further reduce background from photoproduction, y estimated from the positron energy 
and angle is required to be y e < 0.95. 

• The net transverse momentum Pt is expected to be close to zero and is measured with an error 
approximately proportional to \JEt- To remove cosmic rays and beam-related background, 
Pt is required to be less than A\/ Et (GeV). 

The efficiency of these cuts for selecting DIS events with Q 2 > 400 GeV 2 as determined by MC 
is, on average, 80%. It is approximately uniform over the kinematic phase space except for the 
region of high y and low Q 2 where the efficiency decreases due to the positron energy requirement. 
Uncertainties in the simulation of the efficiency arising from the diffractive contribution to the 
cross-section which is not included in the MC are neglected since the diffractive contribution is 
small for Q 2 > 400 GeV 2 . 

6.3 Final event sample 

After the event selection 37438 events with Q 2 A > 400 GeV 2 remain. Distributions from data and 
signal MC are compared in Fig. |||. Shown are the following: the positron energy; the momentum of 
the positron track; the energy not assigned to the positron in an rj — $ cone of 0.8; the Z position 
of the event vertex; the ^-distribution; the scattering angle of the positron, 8 e ; and 7^,, the angle 
of the hadronic system as obtained from (|i~2|). Good agreement between data and MC is seen in 
most variables. Slight disagreements between data and MC at lower positron energies and at low 
values of 7^ may indicate simulation errors in either fragmentation or detector response and are 
accounted for in the systematic uncertainties (see Sect. |8|). 

The photoproduction background is estimated to be less than 0.5% over the full phase space and 
less than 3% at high y. Background from prompt photon events is negligible. 

Backgrounds not related to ep reactions, such as cosmic rays or beam-related background, are inves- 
tigated by studying events in the tails of distributions of the calorimeter timing and of Pt / y/Er, 
and events from unpaired positron or proton bunches. No evidence for such background is ob- 
served and an upper limit of 0.05% at low Q 2 , rising to 0.6% at high Q 2 , is placed on any such 
contamination. 



7 Cross-section determination 

The single differential cross-sections are determined using bin-by-bin unfolding. The measured 
cross-section in a particular bin, <7 meas , is determined from 

Cmeas — ~A£, ' 

where iV b s is the number of observed events in the bin, N^ g is the estimated number of background 
events, A is the acceptance and C is the integrated luminosity. The acceptance, defined from the 
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MC as the number of events reconstructed within the bin divided by the number of events generated 
in that bin, derives from the selection efficiency and the purity in the bin. 

The measured cross-section includes the radiative effects discussed in Sect. The correction factor 
to provide the Born level cross-section is defined as 

o-SM 

°rad 

The numerator is obtained by numerically integrating (|l|) over the bin with fixed a = 1/137, Mi 
= 91.175 GeV, and s'm 2 6 w = 0.232. The value of ofld' the bin cross-section including radiation, is 
calculated using the HERACLES MC. The measured Born level cross-section is then given by 

^Born — ^measC ra( j. (15) 

Finally, the quoted differential cross-section, for example da/dQ 2 , is calculated as 

do_ OBorn ^jrcn /•. ^ 

dQ 2 < rn dQ 2 ■ [ b) 

A similar procedure is used for dajdx and da/dy. In this manner the effects of all the selection 
cuts are corrected (Sect. ||) and the cross-sections are extrapolated to the full kinematic range. In 
particular the MC is used to extrapolate beyond the y-region restricted by the y e cut. 

For the differential cross-section da/dQ 2 the selected events are divided into 20 bins in log 10 Q 2 A . 
The first 12 of these bins divide the Q range of 400 to 3200 GeV 2 into bins of equal width. The 
remaining 8 bins divide the Q 2 range of 3200 to 51200 GeV 2 using bins that are twice as wide. 
For da/dx, the events are divided in log 10 a; into five bins per decade in the range x < 0.4 for 
Q 2 > 400 GeV 2 , x < 0.63 for Q 2 > 2500 GeV 2 and x < 1 for Q 2 > 10000 GeV 2 . These limits 
restrict the analysis to a region of small bin-to-bin migration. To measure the differential cross- 
section da/dy, the selected events are divided into y bins of width 0.05 for Q 2 > 400 and 2500 GeV 2 
and bins twice as large for Q 2 > 10000 GeV 2 . The values at which da/dQ 2 and da/dx are quoted, 
Q 2 and x c , are chosen to be near the logarithmic center of each bin. The cross-section da/dy is 
quoted at the center y c of each bin. 

Typical bin purities are 75%, where the purity quantifies migration effects and is defined as the 
ratio of the number of events generated and measured in the bin to the number of events measured 
in the bin. 

For Q 2 > 400 GeV 2 , the Q 2 A resolution given by the r.m.s. of (Q 2 A — Q 2 )/Q 2 is approximately 
constant at 5%. Q 2 A is slightly biased towards higher Q 2 values, mainly due to the effect of 
initial-state QED radiation. 

The distribution in (x DA — x)/x has an r.m.s. of 9% for 0.03 < x < 0.2. In the remaining part of 
the considered x range the r.m.s. increases to 12%. The distribution in (y DA — y) has an r.m.s. of 
0.05, independent of y. 

The statistical errors are calculated using the square root of the number of measured events N for 
N > 100 and otherwise from 68% Poisson confidence intervals around N. 

The values of Q 2 , x c and y c , the number of observed events, iV bs, the estimated number of pho- 
toproduction background events, iV"b g , the acceptance, A, and the radiative correction factor, C ra d, 
are given m Tables | to g 
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8 Systematic uncertainties 



Systematic uncertainties associated with detector effects are estimated by re-calculating the cross- 
section after re- weighting and modifying the MC to account for discrepancies between data and MC 
in reconstruction and efficiencies. Cut values are varied where this method is not applicable. The 
individual uncertainties are added in quadrature separately for the positive and negative deviations 
from the nominal cross-section values to obtain the total systematic uncertainty. The uncertainty 
on the luminosity of the combined 1994-1997 sample is 1.6% and is not included in the total 
systematic uncertainty. The remaining uncertainties are discussed in detail belowQ: 



The main uncertainty in the trigger chain is expected to come from the first level. Re- 
weighting the MC efficiency to that derived from data results in systematic uncertainties 
below 1%. 

The positron identification efficiency predicted by the MC is checked with a data sample of 
NC DIS events selected using independent requirements such as high Et in the trigger and an 
isolated high-p f track which is associated with the scattered positron. The efficiency curves 
from MC and data agree to better than 0.3% without a bias. For lower positron energies, 
the efficiency is checked using elastic QED Compton scattering events (see Sect. |6.2| ). The 
difference in the efficiencies between data and MC is found to be smaller than 1.5% at the 
smallest electron energy. In addition, loose cuts for forward, high-energy positrons beyond 
the tracking acceptance are used to select candidates which are then inspected. The results 
are consistent with the MC. 

The uncertainty in the positron energy scale (as described in Sect. results in systematic 
variations in the cross-section of 1% at high y and in negligible uncertainties otherwise. Vary- 
ing the positron identification efficiency according to the elastic QED Compton scattering 
events and the track finding efficiency, as derived from a comparison between data and MC, 
produces uncertainties of at most 2%. The positron isolation requirement is estimated by 
varying the cut value by ± 2 GeV and causes systematic uncertainties of up to 2%. Uncer- 
tainties in the measurement of the positron angle produce a systematic variation of up to 2%; 
not using the fiducial cut resulted in variations of less than 2%. 

The uncertainty associated with the reconstruction of 7^ is investigated as follows: by varying 
the calorimeter energy scale separately for RCAL, BCAL, FCAL according to the energy scale 
uncertainty described in Sect. |5.5|; by varying the 7 max parameter in the correction described 



in Sect. 5.4 in a range still compatible with an optimal reconstruction of 7^; by changing the 
energy of the calorimeter cells adjacent to the forward beampipe (and not associated with 
the current jet) based on the uncertainty estimated from a data-MC comparison reflecting 
uncertainties in the simulation of the proton remnant; by excluding events with r jh< 0.1 rad 
to check for the effect of loss of hadronic energy through the forward beam hole of the 



calorimeter; by not using the modified hadronic energy resolution in the MC (see Sect. p.5| 



and by exploring the differences between predictions from the LEPTO and ARIADNE [41] 
models of fragmentation. The last mentioned effect gives the dominating contribution to the 
systematic uncertainty. The net result is an estimated systematic uncertainty of less than 3% 



3 Note that the limits of error referred to are on the absolute magnitude of the error. 
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in the differential cross-sections at low Q 2 and low x, increasing to approximately 8% at high 
Q 2 or high x. 



• The uncertainty arising from the limited knowledge of the shape of the vertex distribution in 
the Z coordinate (see Sect. ||) is at most 1%. 

• Systematic uncertainties due to background removal are estimated by varying the cuts on 
5 and y e in a range that changes the expected background by more than 10% and varying 
the cut on / V Et such that signal events are strongly affected. The uncertainties in the 
cross-section are below 2% for most of the kinematic range; they increase to 8% at high Q 2 
due to the y e cut. The systematic uncertainty arising from a possible underestimation of the 
photoproduction background is obtained from the effect of doubling the background predicted 
from the MC. This results in negligible changes in the cross-sections for most of the kinematic 
range, except at high y where a change of 3% is found. 



An important cross check is the determination of the cross-section using the positron variables [45] 
rather than the double- angle variables. The results from the two methods agree to better than 2% 
for all points. 



9 Results 

The differential cross-sections for NC scattering, da/dQ 2 , da/dx and da/dy are presented in Figs. [| 
to and Tables |l| to || as functions of Q 2 , x and y, respectively. The cross-section da/dQ 2 decreases 
by six orders of magnitude between Q 2 = 400 and 40000 GeV 2 . This decrease is dominated by the 
photon propagator. The cross-section da/dx is shown for different Q regions, Q 2 above 400, 2500 
and 10000 GeV 2 , respectively. A slow fall-off is observed towards x = 0.5 followed by a rapid drop 
towards x = 1. The selection Q 2 > 10000 GeV 2 limits the NC process by kinematics to the region 
x > 0.1 where the contribution from valence quarks (u v ,d v ) is expected to dominate. The cross- 
section da/dy is presented for the same regions in Q 2 as used for da/dx. For Q 2 > 400 GeV 2 the 
bulk of the cross-section is concentrated at small values of y. For Q 2 > 10000 GeV 2 the cross-section 
is approximately constant with y. 

The predictions of the Standard Model (solid curves with PDF uncertainties, see Sect. ||) give 
a good description of all measured cross-sections, except for da/dQ 2 in the highest bin with 
Q 2 > 36200 GeV 2 where two events are observed while 0.27 are predicted by the SM. These events 
were reported previously H as part of an excess seen at high x and high y, obtained from the 
first half of the data. No additional events were observed in the high-Q 2 bin after doubling the 
integrated luminosity. As mentioned earlier, the present analysis is not optimized for the detection 
of narrow high-mass lepton-hadron resonances; an analysis of this type is in progress. 

NC scattering at high Q 2 is sensitive to the contribution from the Z°. According to the SM, the Z° 
contribution reduces the cross- section for Q 2 > 10000 GeV 2 by about 25%. The presence of the Z 
contribution in NC deep inelastic scattering was first demonstrated at Q 2 around 1-2 GeV 2 through 



the observation of an asymmetry in the scattering of polarized electrons on deuterons [47]. The high 
precision data from the present analysis provide sensitivity in the inclusive NC DIS cross-section to 
the Z° contribution at high Q 2 . This is of particular interest as this measurement in the space- like 
region is complementary to the time-like production of the Z° in pp and e + e~ annihilation and thus 
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is an important test of the Standard Model. In Fig. |8] the measured cross-sections are compared with 
the SM predictions, varying the mass of Z° in the propag ator (see ([§)), to values of M z = 40, 91 
GeV and infinity, while keeping the couplings fixed. Figure |8|(a) shows the ratio of the measured 
cross-section da/dQ 2 to the prediction of the SM and Fig. |8](b) presents da/dx as a function of x for 
Q 2 > 10000 GeV 2 . After separately normalizing the SM prediction with or without Z° exchange 
to the data in the full Q 2 range available (and thus essentially eliminating uncertainties arising 
from either the luminosity measurement or the PDFs), x 2 values are calculated from the cross- 
sections in a Q range sensitive to Z°-exchange {da/dQ 2 for 2100 <Q 2 < 10000 GeV 2 and da/dx 
for Q 2 > 10000 GeV 2 ). Considering statistical errors only, the SM prediction yields x 2 = 10-3 for 
10 degrees of freedom, corresponding to a probability prob(SM)= 41%. In contrast, omitting the 
Z° contribution to the cross-section, yields x 2 = 26.4, i.e. prob(SM without Z°)= 0.3%. Taking 
into account each source of systematic uncertainties at a time induces x 2 variations in the range 
8.9-11.8 (22-29) for the SM prediction with (without) Z° exchange, implying prob(SM without 
Z°)< 1.4%. 

10 Summary 

We have studied deep inelastic neutral-current e + p scattering based on data collected during 
1994-1997 with a total luminosity of 47.7 pb _1 . The differential cross-sections da/dQ 2 (for 
400 < Q 2 < 51200 GeV 2 ) and da/dx, da/dy (for Q 2 > 400, 2500, 10000 GeV 2 ) have been mea- 
sured with typical statistical and systematic errors of 3-5% for Q 2 < 10000 GeV 2 . The cross-section 
da/dQ 2 falls by six orders of magnitude between Q 2 = 400 and 40000 GeV 2 . The predictions of 
the Standard Model are in very good agreement with the data. Complementing the observations of 
time-like Z° contributions to fermion-antifermion annihilation, these data provide direct evidence 
for the presence of Z° exchange in the space-like region explored by deep inelastic scattering. 
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Table 1: The differential cross-section da/dQ 2 for the reaction e + p — ► e + X. The following quantities 
are given for each bin: the Q 2 range; the value at which the cross-section is quoted, Q 2 ; the number 
of selected events, iV b s ; the number of expected background events, iVb g ; the acceptance, A; the 
radiative correction factor, C ra( j (see Sect. 0); the measured Born-level cross-section da/dQ 2 ; and 
the Born-level cross-section predicted by the Standard Model using CTEQ4D parton momentum 
distributions. The first error of each measured cross-section gives the statistical error, the second 
the systematic uncertainty. 



Q 2 


range 




Ql 


N ohs 




A 




da/dQ 2 (pb/GeV 2 ) 


(GeV 2 ) 




(GeV 2 ) 


measured 


SM 


400.0 


475 


7 


440 


8504 


2 


4 


0.79 


0.94 


2 


753 ±0.035 


+0.066 
—0.051 




2.673 




475.7 


565 


7 


520 


6319 


2 


4 


0.77 


0.92 


1 


753 ±0.024 


+0.047 
—0.039 




1.775 




565.7 


672 


7 


620 


5008 


2 


8 


0.76 


0.94 


1 


187 ±0.018 


+0.022 
—0.023 




1.149 




672.7 


- 800 





730 


3951 


3 


2 


0.80 


0.94 


(7 


71 


±0.13 


+0.14n 
—0.36/ 


■10- 1 


7.65 ■ 10" 


1 


800.0 


951 


4 


870 


3210 


6 


7 


0.87 


0.93 


(4 


79 


±0.09 


+0 '2l) 


■nr 1 


4.93 ■ 10" 


1 


951.4 


1131 


4 


1040 


2641 


3 


2 


0.89 


0.94 


(3 


21 


±0.07 


+0.06N 
-0.06) 


■10- 1 


3.13 • 10" 


1 


1131.4 


- 1345 


4 


1230 


2000 


1 


2 


0.90 


0.92 


(2 


01 


±0.05 


+0.04\ 
-0.0ZI 


■nr 1 


2.04 ■ 10" 


1 


1345.4 


1600 





1470 


1531 


1 


6 


0.91 


0.93 


(1 


27 


±0.03 


+0.03 \ 
-0.02,! 


■nr 1 


1.28 ■ 10" 


1 


1600.0 


1902 


7 


1740 


1204 


2 





0.91 


0.92 


(8 


49 


±0.26 


+0.17S 
-0.30-1 


■10~ 2 


8.26 ■ 10" 


2 


1902.7 


- 2262 


8 


2100 


863 





4 


0.91 


0.93 


(4 


97 


±0.18 


+ 0.11N 

-0.16- 1 


■10" 2 


5.01 ■ 10" 


2 


2262.8 


- 2690 


9 


2500 


629 





4 


0.92 


0.94 


(3 


05 


±0.13 


+0.06\ 
-0.14/ 


■10~ 2 


3.13 ■ 10" 


2 


2690.9 


- 3200 





2900 


455 


1 


2 


0.91 


0.94 


(1 


99 


±0.10 


+0.07n 
-0.09/ 1 


■10" 2 


2.09 • 10" 


2 


3200.0 


- 4525 


5 


3800 


565 


3 


1 


0.91 


0.93 


(9 


00 


±0.39 


+0.20\ 
-0.241 


■10" 3 


9.77 ■ 10" 


3 


4525.5 


- 6400 





5400 


303 








0.91 


0.91 


(3 


30 


±0.19 


+0.17\ 

-0.10) 


■10" 3 


3.49 ■ 10" 


3 


6400.0 


- 9050 





7600 


162 








0.90 


0.93 


(1 


32 


±0.10 


+0.02\ 
-0.07) 


■10" 3 


1.20 ■ 10" 


3 


9050.0 


- 12800 





10800 


63 








0.86 


0.93 


(3 


69 


+0.53 
-0.47 


+0.08\ 
-0.11) 


■10" 4 


3.64 ■ 10" 


4 


12800.0 


- 18102 





15200 


20 








0.81 


0.93 


(8 


9 


+2.5 
-2.0 


+0.7\ 
-0.6) 


■10" 5 


10.0 ■ 10" 


5 


18102.0 


- 25600 





21500 


8 








0.86 


0.96 


(2 


4 


+ 1.2 
-0.8 


+0.4\ 

-0.1J 


■10" 5 


2.2 ■ 10" 5 




25600.0 


- 36203 





30400 











0.86 


0.90 






< 6.0 ■ 


10" 6 




3.7- 10" 6 




36203.0 


- 51200 





43100 


2 








0.89 


1.00 


(2 


6 


+3.5 
-1.7 


+0.7\ 
-0.2) 


•10" 6 


0.4 ■ 10" 6 
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Table 2: The differential cross-section da/dx for the reaction e + p — ► e + X for Q 2 > 400, 2500, and 
10000 GeV 2 . The following quantities are given for each bin: the Q 2 and x range; the value at 
which the cross-section is quoted, x c ; the number of selected events, iV b s ; the number of expected 
background events, A^b g ; the acceptance, A; the radiative correction factor, C ra d (see Sect. 0); 
the measured Born-level cross-section da/dx; and the Born-level cross-section predicted by the 
Standard Model using CTEQ4D parton momentum distributions. The first error of each measured 
cross-section gives the statistical error, the second the systematic uncertainty. 



o 2 ■ 

mm 


X 


range 






N obs 




A 




da/dx (pb) 


(GeV 2 ) 




measured 


SM 


400 


(0.63 - 


1.00) -10" 


2 


0.00794 


2307 


6.3 


0.67 


0.98 


(1.96 ±0.05 ±0.05) 


10 4 


1.88 


■ 10 4 




(0.10 - 


0.16) -10" 


1 


0.0126 


4352 


11.0 


0.84 


0.94 


(1.78 ±0.03 t° ° l) 


10 4 


1.70 


• 10 4 




(0.16 - 


0.25) -10" 


1 


0.0200 


5026 


6.7 


0.85 


0.94 


(1.27 ±0.02 t°m) 


10 4 


1.24 


• 10 4 




(0.25 - 


0.40) -10" 


1 


0.0316 


5283 


2.7 


0.87 


0.92 


(8.08 ±0.12 +°iq) 


10 3 


8.04 


• 10 3 




(0.40 - 


0.63) -10" 


1 


0.0501 


5028 


2.7 


0.87 


0.92 


(4.83 ±0.07 t :°l) 


10 3 


4.95 


• 10 3 




(0.63 - 


1.00) -10" 


1 


0.0794 


4782 


0.0 


0.86 


0.93 


(2.96 ±0.05 ±Sloe) 


10 3 


2.92 


• 10 3 




0.10 - 


0.16 




0.126 


4219 


0.4 


0.86 


0.92 


(1.63 ±0.03 torn) 


10 3 


1.66 


■ 10 3 




0.16 - 


0.25 




0.200 


3512 


0.0 


0.82 


0.92 


(9.04 ±0.17 tg.lr) 


10 2 


8.68 


■ 10 2 




0.25 - 


0.40 




0.316 


2276 


0.4 


0.87 


0.92 


(3.51 ±0.08 to.w) 


10 2 


3.72 


■ 10 2 


2500 


(0.25 - 


0.40) 10" 


1 


0.0316 


58 


1.2 


0.72 


1.01 


f 1 1 c +0.18 +0.09\ 
yi.LO _ 16 _o.l3/ 


10 2 


1.19 


• 10 2 




(0.40 - 


0.63) -10" 


1 


0.0501 


252 


2.7 


0.91 


0.95 


(2.40 ±0.16 to w) 


10 2 


2.49 


• 10 2 




(0.63 - 


1.00) -10" 


1 


0.0794 


340 


0.0 


0.94 


0.93 


(1.94 ±0.11±°;°I) 


10 2 


2.16 


• 10 2 




0.10 - 


0.16 




0.126 


421 


0.0 


0.93 


0.91 


(1.51 ±0.08 t° °ol) 


10 2 


1.51 


• 10 2 




0.16 - 


0.25 




0.200 


356 


0.0 


0.93 


0.92 


(8.16 ±0.44 to.lt) 


10 1 


8.97 


■10 1 




0.25 - 


0.40 




0.316 


265 


0.4 


0.87 


0.91 


(4.02 ±0.25 ±0.20) 


10 1 


4.15 


■10 1 




0.40 - 


0.63 




0.501 


112 


0.0 


0.84 


0.93 


(1.10 ±0.11 jftE) 


10 1 


1.06 


■10 1 


10000 


0.10 - 


0.16 




0.126 


7 


0.0 


0.62 


1.06 


A r +2.4 +1.0 
-1.7 -0.6 




3.3 






0.16 - 


0.25 




0.200 


19 


0.0 


0.82 


0.94 


r n +1.4 +0.4 
d - u -1.2 -0.5 




6.0 






0.25 - 


0.40 




0.316 


23 


0.0 


0.88 


0.91 


c +0.9 +0.1 
°-° -0.7 -0.2 




4.1 






0.40 - 


0.63 




0.501 


12 


0.0 


0.85 


0.93 


1 +0.5 +0.2 
l ' £ -0.3 -0.1 




1.3 






0.63 - 


1.00 




0.794 


2 


0.0 


0.92 


0.95 


(5-4 tH til) 


10" 2 


3.6 ■ 


10~ 2 
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Table 3: The differential cross-section da/dy for the reaction e + p — > e + X for Q 2 > 400 GeV 2 . The 
following quantities are given for each bin: the Q 2 and y range; the value at which the cross-section 
is quoted, y c \ the number of selected events, iV b s ; the number of expected background events, N^ g ; 
the acceptance, A; the radiative correction factor, C ra d (see Sect. |7|); the measured Born-level cross- 
section da/dy; and the Born-level cross-section predicted by the Standard Model using CTEQ4D 
parton momentum distributions. The first error of each measured cross-section gives the statistical 
error, the second the systematic uncertainty. 



Qmin 


y range 


Vc 


N ohs 




A 


(--rad 


da/dy (pb) 


(GeV 2 ) 


measured 


SM 


400 


0.00 


- 0.05 


0.025 


5613 


0.0 


0.80 


0.95 


(3.82 ±0.06 t°'lo) 


•10 3 


3.87 


10 3 




0.05 


- 0.10 


0.075 


5844 


0.4 


0.84 


0.94 


(2.68 ±0.04 t°'it) 


■10 3 


2.64 


10 3 




0.10 


- 0.15 


0.125 


4128 


0.4 


0.86 


0.93 


(1.86 ±0.03 ±004) 


■10 3 


1.89 


10 3 




0.15 


- 0.20 


0.175 


3231 


0.0 


0.88 


0.93 


(1.43 ±0.03 i^oi) 


■10 3 


1.47 


10 3 




0.20 


- 0.25 


0.225 


2685 


0.0 


0.86 


0.93 


(1.22 ±0.03 


■10 3 


1.20 


10 3 




0.25 


- 0.30 


0.275 


2226 


0.0 


0.85 


0.92 


(1.01 ±0.02+^4) 


■10 3 


1.01 


10 3 




0.30 


- 0.35 


0.325 


1939 


0.0 


0.85 


0.92 


(8.83 ±0.22 to™) 


■10 2 


8.68 


10 2 




0.35 


- 0.40 


0.375 


1731 


2.4 


0.87 


0.93 


(7.76 ±0.21 to.lt) 


■10 2 


7.57 


10 2 




0.40 


- 0.45 


0.425 


1547 


1.6 


0.86 


0.91 


(6.90±0.19±g;i|) 


■10 2 


6.69 


10 2 




0.45 


- 0.50 


0.475 


1389 


0.0 


0.90 


0.92 


(5.95 ±0.18 toll) 


■10 2 


5.96 


10 2 




0.50 


- 0.55 


0.525 


1308 


1.2 


0.91 


0.91 


(5.44 ±0.17 1°-™) 


■10 2 


5.36 


10 2 




0.55 


- 0.60 


0.575 


1205 


3.6 


0.91 


0.91 


(5.00 ±0.16 ±^7) 


■10 2 


4.85 


10 2 




0.60 


- 0.65 


0.625 


1104 


0.4 


0.92 


0.92 


(4.64 ±0.15 to.m) 


■10 2 


4.43 


10 2 




0.65 


- 0.70 


0.675 


978 


0.8 


0.89 


0.91 


(4.20 ±0.15 t o.l 3 ) 


■10 2 


4.06 


10 2 




0.70 


- 0.75 


0.725 


849 


2.8 


0.84 


0.94 


(3.99 ±0.15 t°oi° 5 ) 


■10 2 


3.75 


10 2 




0.75 


- 0.80 


0.775 


676 


1.2 


0.74 


0.92 


(3.52 ±0.15 to.H) 


■10 2 


3.48 


10 2 




0.80 


- 0.85 


0.825 


460 


2.4 


0.55 


0.93 


(3.28 ±0.16 to.lt) 


■10 2 


3.25 


10 2 




0.85 


- 0.90 


0.875 


311 


5.1 


0.41 


0.94 


(2.98 ±0.18 toil) 


■10 2 


3.06 


10 2 




0.90 


- 0.95 


0.925 


209 


5.9 


0.24 


0.96 


(3.35 ±0.25 toll) 


■10 2 


2.90 


10 2 
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Table 4: Continuation of Table pc The tabulation of results for the differential cross-section da/dy 
for Q 2 > 2500 and 10000 GeV 2 . 



o 2 ■ 

mm 


y range 
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AT 
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A. 


n 

t-rad 


da/dy (pb) 


(GeV 2 ) 
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SM 
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0.05 
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+0.13\ 
-0.50) 


■10 1 
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0.95 


(7.55 ±0.59 


+0.23N 
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6.58 
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0.425 
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4.70 
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(A 04 + ' 46 
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■10 1 
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85 
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0.89 
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3.82 
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0.55 


- 0.60 
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86 
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/o 4 i +0.42 


+0.08\ 
-0.17-' 
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3.48 
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- 0.65 
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72 
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0.90 


(n Q9 +0.40 
_ .35 


+0.15\ 
-0.06^ 
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3.18 


■10 1 




0.65 


- 0.70 


0.675 


72 
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0.92 


0.89 


(2 94 + - 40 
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+0.13\ 
-0.15/ 


■10 1 


2.93 


■10 1 




0.70 


- 0.75 


0.725 


62 
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0.94 


0.91 


O cq +0.37 
l i - tJJ -0.33 


+0.27\ 
-0.06^ 


■10 1 


2.72 


■10 1 




0.75 


- 0.80 


0.775 


74 


0.0 


0.93 


0.92 


(3.08l° ;« 


+0.08\ 
-0.21/ 


■10 1 


2.54 


■10 1 




0.80 


- 0.85 


0.825 


44 


1.6 


0.85 


0.93 


(i.93lg:IS 


+0.14n 
-0.12 1 


■10 1 


2.39 


■10 1 




0.85 


- 0.90 


0.875 


44 


1.2 


0.83 


0.92 


(1 01 + 037 

r- ul -0.32 


+0.26\ 
-0.il) 


■10 1 


2.27 


■10 1 




0.90 


- 0.95 


0.925 


55 


1.2 


0.74 


0.93 


(0 Se; +0.46 
-0.41 


+0.16\ 
-0.23/ 


■10 1 


2.17 


■10 1 


10000 
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1.08 


(>7 o +18.0 
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•10" 2 


5.0 • 
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1 - Z -0.5 


+0.1 

-0.2 




1.3 






0.3 


- 0.4 


0.35 
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0.0 


0.82 
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-0.3 




2.6 
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0.55 
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+0.0 
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2.4 






0.6 


- 0.7 


0.65 


10 


0.0 


0.86 


0.91 


o 9 +1.0 
-0.7 


+0.2 
-0.0 




2.2 
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- 0.8 
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10 


0.0 


0.95 


0.92 


2 +0 ' 9 

z - u -0.6 


+0.2 

-0.1 




2.0 
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- 0.9 


0.85 


9 


0.0 


0.93 


0.95 


1 q +0-9 
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-0.1 




1.8 






0.9 
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0.95 
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0.31 
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2.8 til 
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Figure 1: The effect of the correction used for the reconstruction of the angle 7^ as described in 
Sect. |5.4| . (a) and (b): the ratio of total energy E rem removed by the correction and the hadronic 
transverse momentum Pt,h divided by sin 7/1 for two different ranges in 7^. The improvement in 
the reconstruction of 7^ by applying the correction is shown for the relative bias in 7^ (c) and the 
relative resolution in 7^ (d), where relative refers to the normalization to the true 7; both are shown 



as a function of the true y as defined in (11). The result with (without) correction as obtained from 
a NC MC is shown as filled (unfilled) symbols. 
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Figure 2: (a) the ratio of the positron energy E' e (r&w) measured in BCAL to that calculated from 
double-angle variables, E^ A , as a function of the Z-position of the positron in BCAL. The dips 
coincide with the cell boundaries in BCAL. (b) the same after applying a non-uniformity correction 
as described in Sect. 5.5, yielding Eg (cor). Open circles show MC and dots show data. 
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Figure 3: A comparison of data (dots) and MC (histogram) for the ratio of the transverse momen- 
tum of the hadronic final state, Pr,h, and the positron, Pr, e , for Pt,h > 20 GeV and 7^ pointing 
into the FCAL (a) and into the BCAL (b). 
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ZEUS NC 1994 - 97 
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Figure 4: A comparison of data (points) and the signal MC predictions (histograms) for 

(a) E' e , the energy of the scattered positron; the inset shows the high energy part of the distribution; 

(b) Pg rk , the momentum of the positron track; (c) the energy E cone not assigned to the positron 
in an r]-& cone of 0.8; (d) the Z position of the event vertex; (e) 5 = T,(Ei — p z i)', (f) &e, the angle 
of the positron; and (g) j^, the angle of the hadronic system. Only events which pass all event 
selection cuts are plotted in (e) and (f); in (a)-(e), events are plotted which pass all other event 
selection cuts except those on the variable displayed; the cuts to be applied to this variable are 
indicated by the vertical lines. 
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Figure 5: The high-Q 2 e + p NC DIS cross-section, da/dQ 2 , for data (points with error bars) and 
the Standard Model predictions using the CTEQ4D parton momentum distributions (a). Also 
plotted is the ratio of data to the prediction (b). The inner error bars (delimited by the horizontal 
lines) show the statistical errors, the outer ones the statistical and systematic uncertainties added 
in quadrature. The shaded region gives the uncertainty in the Standard Model prediction due to 
the uncertainty in the parton momentum distributions (PDF) |2(J . 
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Figure 6: The high-Q 2 e + p NC DIS cross-section, da/dx, for data (points with error bars) and 
the Standard Model predictions using the CTEQ4D parton momentum distributions. Plotted 
are cross-sections for (a) Q 2 > 400 GeV 2 , (b) Q 2 > 2500 GeV 2 , and (c) Q 2 > 10000 GeV 2 . The 
inner error bars (delimited by the horizontal lines) show the statistical errors, the outer ones 
the statistical and systematic uncertainties added in quadrature. The shaded region gives the 
uncertainty in the Standard Model prediction due to the uncertainty in the parton momentum 
distributions (PDF) [p0|. The insets show the ratio of data to the prediction. 
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Figure 7: The high-Q 2 e + p NC DIS cross-section, da/dy, for data (points with error bars) and 
the Standard Model predictions using the CTEQ4D parton momentum distributions. Plotted 
are cross-sections for (a) Q 2 > 400 GeV 2 , (b) Q 2 > 2500 GeV 2 , and (c) Q 2 > 10000 GeV 2 . The 
inner error bars (delimited by the horizontal lines) show the statistical errors, the outer ones 
the statistical and systematic uncertainties added in quadrature. The shaded region gives the 
uncertainty in the Standard Model prediction due to the uncertainty in the parton momentum 
distributions (PDF) ^0|. The inset in (a) shows the ratio of data to the prediction. 



27 



ZEUS NC 1994 - 97 



O 

ns 10 



^ 10 
3 9 

£ 8 

^ 7 
* 6 
5 
4 
3 
2 
1 



1.2 



0.8 




a) 



10 



10 



10 



10 



Q 2 (GeV 2 ) 









• ZEUS 94-97 e + p NC 




* 

* 


** %- 

x 

X 

N 

\ 

s 


M z = 00 






s 

i. V 

\ 

i \ 


— M z = 91GeV 


= < 


' s < 


» 


M z = 40 GeV 


i 

Z ' 

i j 








- ' /' 

- '/ 

- i J 

- '/ 

- if 


Q 2 >10000 GeV 2 


b) 








1 1 1 ^ A 1 



0.1 



0.2 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 



Figure 8: The points with error bars give the ratio of the measured cross-sections, da/dQ 2 , to the 
SM prediction using the CTEQ4D parton momentum distributions and fixing the Z° mass, Mz, 
at its nominal value of 91.175 GeV (a) and da/dx for Q 2 > 10000 GeV 2 (b). The three lines show 
the SM predictions for M z = 91.175 GeV (solid line), for M z = 40 GeV (dotted line) and for no 
Z° contribution (dashed line). The inner error bars (delimited by the horizontal lines) show the 
statistical errors, the outer ones the statistical and systematic uncertainties added in quadrature. 
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